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A summary of the extraction of light quark masses from both QCD sumrules and lat- 
tice QCD simulations is presented. The focus is on providing a careful statement of 
the potential weaknesses in each calculation, and on integrating the work of different 
collaborations to provide a coherent picture. 

1. Introduction 

Significant progress has been made since 1996 in the determination of the light 
quark masses from both lattice QCD (LQCD) and QCD sumrules (SR). The evolu- 
tion in time of the ranges quoted by the Particle Data GroupEJ and our best estimates 
as of DPF2000 are given below. Since the scale dependence between 2 and 1 GeV 
is very large, we quote all results in the MS scheme at scale /x = 2 GeV. Results 
at the lower scale 1 GeV may be obtained via mq{\ GeV) = 1.38 mq(2 GeV) where 
the factor 1.38 corresponds to 4-loop running, with as{m'^) — 0.3340. In quoting 
separate values for m„ and vrid we have employed the ratio mu/md = 0.553, ob- 
tained from ChPT. The basis for these estimates is outlined briefly bdpw. A more 
detailed discussion will appear in an extended version of this writeup □. 

1996 2000 Proposed SR Proposed LQCD 

niu 2-8 MeV 1 - 5 MeV 2.4 - 3.8 MeV 2.2 - 2.7 MeV 

rud 5-15 MeV 3-9 MeV 4.3 - 6.9 MeV 3.8 - 4.9 MeV 

ms 100 -300 MeV 75 - 170 MeV 83 - 130 MeV 78 - 100 MeV 
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2. QCD Sum Rule Determinations 

We focus on recent (>1995) extractions of m„ + and mg, which employ either 
Borel (Laplace) transformed sum rules (BSR's) or finite energy sum rules (FESR's). 
For a typical correlator 11 with spectral function p, BSR's approximate p(s) by 
its OPE version for s > sq. Uncertainties in the the corresponding "continuum" 
spectral contribution are small only if sq/M'^ is significantly > 1 for Borel masses, 
M, in the "stability window" of the analysis. In contrast, the spectral ansatz has an 
unphysical gap if sq lies significantly beyond the squared mass of the last resonance 
appearing in the low-s part of the ansatz. FESR's relate the integral from to so 
of p{s)w{s) (with w{s) any analytic function) to the integral of n(s)i(;(s) over the 
circle \s\ — sq in the complex s-plane. Using the OPE representation for 11 in the 
latter integral is a potential problem since the OPE is expected to break down near 
the real timelike axis. A study of the isovector vector (IVV) channel using hadronic 
r decay dataQ shows that FESR's based on weights w{s) — s*^ are not well-satisfied 
at scales ~ 2 — 3 GeV^ but those based on "pinched weights" (those satisfying 
w(so) = 0), which we will denote PFESR's, are, in contrast, very well-satisfied. 

2.1. The Isovector and Isospinor Pseudoscalar Sum Rules 

Since = (m^ -I- rrij) : qii^zqj'-, rriu + ma and -t- can be determined 

from sjim rules for the corrissponding ij — ud andrfis correlators. Recent analyses 
(BPREi and its update, P98EI, for ij = ud and DPSQ for ij = us) employ the 4-loop 
]D — OPE expression. Known — 4,6 contributions are at the few percent level, 
while other non-perturbative effects are assumed to be small. The tt (K) spectral 
contributions are known, but those from the (excited) resonance region are not. A 
modified sum of two Breit-Wigners is used to model the latter, the overall scale 
being set by normalizing the sum of resonance tails at threshold to values given by 
ChPT. P98 uses non-PFESR's (w(s) = 1 and s) and so ~ 2 ^ 3.5 GeV^; DPS use 
the BSR framework and, for A-qcjj = 380 MeV, so ~ 6 ^ 8 GeV^. The resuhs are 

[m.„ + mrf](2 GeV) = 9.8 ± 1.9 Mevi and m,.(2 GeV) = 112 ± 18 Mevi. PFESR's 
employing the BPR/P98 ansatze (tuned originally using non-PFESR 's) in the same 

range sch~ 2-3 GeV^ used by BPR/P98 yield very poor OPE/spectral integral 
matched. Possible sources of this problem are the resonance normalization pre- 
scription and the use of non-PFESR's at scales where they are not well-satisfied in 
the IVV channel. Potential problems for DPS are the use of an input assumption 
about the relative strengths of the K{1A60) and i4r(1830) contributions and the 
existence of a spectral gap between s ^ 4 and 6 GeV'^. To shift focus to spectral 
normalization at the resonance peaks, rather than at threshold, an incoherent sum 
of Breit-Wigners can be employecL and a PFESR analysis used to fit both rm + ruj 
and the resonance decay constanta3. The resulting optimized OPE/spectral integral 
matches are at the ~ 1% level or better in the fitting window, 2.8 — > 3.6 GeV^, and 
correspond to [m„ + mrf](2 GeV) = 9.9 ± 1.0 MeV, ms{2 GeV) = 116 ± 5 MeV. 
The above results neglect direct instantons, and other non-perturbative effects not 
present in the OPE. Two examples show that non-trivial uncertainties are associ- 
ated with this neglect. First, a phenomenologically-determined effective tachyonic 
gluon mass-squared, meant to represent additional short-distance non-perturbative 
contributions, lowars the P98 value by 5.6%eI. Second, including instanton liq- 
uid model (ILM)ll3 estimates of direct instanton effects, the PFESR extraction, 
[niu + mrf](2 GeV) = 9.9 MeV is reduced to 7.8 MeVH. 

2.2. The Light- Strange Scalar Sum Rule 

Since = i{ms — niu) :su:, mg — niu — mg can be determined from sum rules 
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for the strange scalar channel. The low-s part of p{s) in this case can be deter- 
mined indirectly, with certain additional theoretical assurnMionSj^om therQmnes 
representation of the timelike Kit scalar formfactor, d^^Ed. JMU, CDPSli3, and 
CPSI13 employed this construction only to fix dKrr at threshold; the tail of a sum- 
of-Breit-Wigners was then normalized to this value, as in the BPR/P98 analyses. 
A purely resonant phase for dxTv, however, yields a slope, d'i^^{s = 0), incompatible 
with ChPTllJ and a poor OPE/spectral integral matchQ; both problems are removed 
if one includes background ccH^pbutions to the phase near tiwcshold-and uses the 
Omnes representation for all slifl'Ej. More recent analyses ( J98li3, M99li3) employ the 
CFNP version of p{s). CFNP and J98 are BSR analyses differing only in the value of 
So, while M99 is a PFESR analysis. CFNP has a good stability plateau for nis, but 
at the cost of a gap in p. J98 has no spectral gap, but also no stability plateau, and 
also potentially non-trivial continuum contributions. The M99 PFESR's show an 
excellent OPE/spectral integral match. The results, ms{2 GeV) = 91 — > 116 MeV 
for CFNP, 116 ± 22 MeV for J98, and 115 ± 8 MeV for M99, are all compatible. 
Additional errors associated with the assumptions made in writing down the Omnes 
representation are at least '^10 MeV in magnitude. 

2.3. Sum Rules Based on Electromagnetic Hadroproduction Data 

Sum rules based on electromagnetic (EM) hadroproduction data000'0 have-sigj 
nificant uncertainties associated with either^l) isospin- breaking correctionalallj 
(for the flavor 33-88 vecixir current sum rulal3) or (2) deviations from ideal mix- 
ing/Zweig rule violations^ (for the flavor 33-ss vector current sum rulec3). We, 
therefore, consider these sumrules uncompetitive with those based on hadronic t 
decay and defer further discussion to the extended version of this write-up. 

2.4. Sum Rules Based on Hadronic Tau Decay Data 

The ratio of the hadronic t decay rate through the flavor f — ij = ud, us vector 

(V) or axial vector (A) current to the r — > VT-ev^ rate, i?]^^^'*"*, can be written as a 
sum of appropriately-weighted integrals of the J = 0, 1 parts of the corresponding 
hadronic spectral function. Since the ij = ud and us correlators are identical in 
the SU{3)f limit, rescaled differences such as i?r;M(i/|KidP — Rt;iis/\Vus\^ = Ai?T- 
(with Vij the ij CKM matrix element and i?T;ij = R^'^-' + R^'^-') vanish in this 
limit. Experimental data allows access to the ud, us spectral distributions, and 
hence to integrals of such correlator differences, whose OPE representations, for 
large enough sq, should be dominated by the D — 2 term, proportional (neglecting 
m^ d) to TO^. The kinematic weights appearing on both sides of the Ai?,- sum 

rule may be supplemented by a factor (1 — s / rn^Y [s / rn^Y without necessitating a 
J — Ql J = 1 experimental decomposition; the analysis is then said to employ the 
"(A;,/) spectral weight". 

The D = 2 OPE representations of the J = + 1 and J — parts of the ud-us 
correlator are known to 3-loop {0{a'^)) and 4-loop (O(ai^) order, respectively. The 
ud spectral distributions are very accurately determinedH while the V + A sum for 
the us case is-known with errors of 6 — 8% in the K* region and ~ 20 — 30% 
abovcr^he K*M. AUrMalyses employ the ALEPH ud and us data. The CDH98B, 
MT98EJ, and CKP98ca analyses are based on the pEpliminary result, Rr.us — 0.155± 
.006; ALEPHQSa, PP99y, KKEOOEI and KMOOEJ on the 1999 publislied version, 
0.1610 ± .0066B; and DHPPC00C3 on a recent update, 0.1630 ± .OOSTES. Increasing 
Rr,us decreases Rr-ud, which is obtained via R^^ud = Rr-had - Rt-us = [(1 - -Be - 
B^)/Be] — Rt-us- Small increases in B^ and B^El have also lowered Rr-had since the 
earlier ALEPH analyses. The high degree of cancellation between ud and us spectral 
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integrals (typically to better than 10%) means the result for nis is very sensitive to 
exactly which values of Rr-ud, Rt-us have been used. For the same reason, iris is 
also sensitive to small variations in the input values of Jk and \Vus\ (|KisP, e.g., 
differs by 2.6% depending on whether one uses the value from K^^ or that from 
3- family CKM unitarity, combined with \Vud\ as extracted from 0+ — ^ 0+ nuclear 
decays; the difference is non- negligible on the scale of the < 10% ud-us residual). 
Different prescriptions for truncating the D — 2 OPE series and estimating the 
truncation error also are present in the literature. In Table 1 we display the impact 
of converting existing results so as to correspond, where possible, to (1) the common 
truncation and error estimate scheme employed by DHPPCOO and (2) Rt-^us — .163. 
A subtlety is involved in the second of these conversions for non-inclusive analyses 
and/or those employing sq ^ rn^; see the extended version of this write-up for 
details. Since the 0{[a'g-^^'^) estimate in the KKPOO framework (which employs an 
effective running coupling and running masses) is not available to us, the KKPOO 
results have been converted only to reflect the new values of the ud and us integrals. 
Two converted versions are given in each case, one corresponding to the PDG2000 
unitarity-constrained CKM set, = .9479, |Ks| = -2225 (CKMU), the other to 
the PDG2000 non-uuitarity-consteiined set, |Kd| = -9735, |Ks| = .2196 (CKMN). 
Results from CDH98E^ and MT98EJ have been omitted since the former employed an 
older, erroneous expression for the integrated Z? = 2, J = OPE contribution and 
the latter an assumption about thei relation between the ud and us J = spectral 
integrals now known to be invalidEil. Also not quoted are the PP99 and ALEPH99 
results, which have been superceded by the joint DHPPCOO update. The first error 
quoted is experimental, the second theoretical. The reader should bear in mind that 
the integrated D = 2 OPE^epiysJor the J = contribution to the ud-us r decay sum 
rules converges very badlyEaOEa (for sq = the ALEPH99 version, for example, 
is ~ 1 + 0.78 + 0.78-f 0.90 + - • E3). This creates potentially large theoretical errors for 
those ("inclusive") analyses which retain both the J = O-f-l and J — contributions 
(ALEPH99, CKP98, PP99, KKPOO and DHPPCOO). While the identification, and 
hence subtraction, of the tt and K pole J = contributions is unambiguous, at 
present noJ = 0/J=l separation exists for the experimental data in the excited 
resonance region. Analyses which subtract the J = contributions (ALEPH99, 
KMOO), and work with the theoretically much better behaved J = 0-t-l contribution 
(the "0-1-1" approach), thus have instead J — subtraction uncertainties. 

Our summary of estimates based on hadronic r decay is given in Table m along 
with the analysis type. Apart from KMOO, these results neglect possible 1? > 6 
contributions. Because the polynomial coefficients grow rapidly with k for the 
(fc, 0) spectral weights, this neglect becomes less safe as k increases. An extended 
discussion of the issues alluded to above, as well as details on the conversion of the 
individual results to the values quoted, will be given in the longer writeup a. 

2.5. Summary 

We consider results for rris based on sum rules involving hadronic r decay data 
both most reliable and to have the greatest potential for improvement (through im- 
proved us spectral data which should be provided by the B factories). Results from 
other analyses are compatible, but with less controlled uncertainties. The agreement 
between all approaches, however, suggests that theoretical uncertainties are under 
reasonable control. Experimental errors in the different r analyses are necessarily 
strongly correlated, and hence cannot be averaged. Consistency of the different 
versions does, however, allow us to average the central values, and quote as our 
error the smallest of the errors shown (which corresponds to the least strong ud-us 
cancellation) . With this prescription, we obtain for our final sum rule determination 
ms{2 GeV) = 114±16 MeV for the CKMU input set and to^(2 GeV) = 101±18 MeV 
for the CKMN input set. As mentioned above, theoretical uncertainties in m„ -fm^ 
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Table 1. Impact of conversion to common input for the hadronic t decay extractions of ms- 
Results for ms(2 GeV) are in MeV, and correspond to Be = .1783, B^, = .1737, Rt;us = .163. All 
except KKPOO also correspond to the DHPPCOO truncation prescription, and to as(m^) = .334. 
Experimental and theoretical errors have been combined in quadrature. 



ms{2 GeV) (MeV) 


Reference 


Analysis Type 


Original 


CKMU input 


CKMN input 


CKP98 


inclusive, (0, 0) 


145 ± 36 


116 ±31 


99 ±34 


KKPOO 


inclusive, (Oc-9) 


125 ± 28 


120 ± 28 


106 ± 32 


KMOO 


(0 + 1), W2M 


115 ± 17 


110 ± 16 


100 ± 18 


DHPPCOO 


inclusive, (0, 0) 


126 ± 31 


124 ± 32 


106 ± 37 


DHPPCOO 


inclusive, (1, 0) 


115 ± 19 


113 ± 19 


102 ± 21 


DHPPCOO 


inclusive, (2, 0) 


100 ± 21 


99 ± 21 


91 ± 21 



are larger. Using the ChPT ratio 2TOs/[m„ + m^] = 24.4 ± 1.5, the CKMU and 
CKMN TOs results correspond to [m„ + md](2 GeV) = 9.3 ± 1.4 and 8.3 ± 1.6 MeV, 
respectively. For comparison, the P98 estimate is either 9.8 MeV or, if one uses the 
P98 1 GeV value together with our scaling factor, 12.8/1.38 = 9.3 MeV, while the 
PFESR/ILM estimate is 7.8 MeV. For our final summary we choose the conserva- 
tive combined ranges 6.7 — 10.7 MeV for to„ + md and 83 — 130 MeV for m^. These, 
as we show next, are in excellent accord with LQCD estimates. 

3. Lattice QCD: Quenched Results 

LQCD simulations provide estimates for two quantities m = (m„ + rnd)/2, and 
wis. mu and are not obtained separately since current simulations do not include 
EM effects. Pjet|ails of how quark masses are extracted from lattice calculations 
are given in and will not be repeated here. We start with a summary of 

the state-of-the-art results given in Table □ and discuss current shortcomings and 
future prospects. This table also provides information on (i) the lattice action used 
in the simulations, (ii) whether extrapolation to the continuum limit was done, and 
(iii) the states used to fix the quark masses and the lattice spacing a. The errors 
quoted include statistical and those due to renormalization constants and chiral and 
continuum extrapolations. Comparing data from different collaborations suggests 
that the quoted statistical errors are realistic and much larger than finite volume 
corrections. The latter three systematic errors will be discussed later. 

At first sight, these estimates suggest a lack of consistency between different 
lattice calculations. We will attempt to argue that this is not so. The key reason 
is that it is only for QCD with three dynamical light flavors (the physical theory) 
that we expect all simulations to give the same results once the extrapolation to 
the continuum limit (a = to remove discretization errors) has been carried out. 
So, before analysing the relative merits of the different calculations one has to 
address the question of concern to all: how valid is the quenched approximation, 
especially since the resulting theory is not even hermitian? From a practical point 
of view quenching introduces two main limitations. First, the spectrum of the 
quenched theory will notrGoiacide with the experimentally observed one. The second 
limitation, discussed in cj'Cj, and known as the problem of quenched chiral logs, 
is more subtle. In a nutshell the problem is that the quenched ?/ persists as a 
Goldstone boson in the chiral limit and its propagator has a single and double pole, 
consequently the chiral expansion of pseudoscalar meson masses, decay constants, 
and Q**ft*k condensates develop enhanced logarithms that are singular in the chiral 
limit llj'cll. The first limitation implies that there will be no consistent set of quark 
masses which reproduces the observed spectrum in the quenched approximation. 
The second is relevant when extrapolating data to the physical u and d quark 
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Table 2. Recent results for quark masses. 0(a) SW stands for non-perturbative 0(a) improved 
Sheikhholeslami-Wohlert (SW) fermion action, Iwasaki for an improved gauge action, and DWF 
for Domain Wall Fermions. In the last column we give the quantity used to set the lattice spacing 
a and (a — > 0) indicates that results were extrapolated to the continuum limit, ro is defined in 
terms of the force at 0.5 fermi and a(ro) ~ a{Mp). 





Action 


m 


ms{MK) 




scale 1/a 
(GeV) 






3.8(1)(3) 


99(3)(8) 


111(7)(20) 


Mp 


(1997)EI 










(a^O) 


(1998)El 


0(a) SW 


4.5(4) 


111(12) 




Mk* 










w 2.7 


n 

(1999)EI 


0(a) SW 


4.8(5) 


111(9) 




Mk- 










« 2.7 


JLQCp. 


Staggered 


4.2o(29 ) 


1U0(/) 


1 on 1 o^ 
129(12 ) 


Mp 
(a 0) 


(1999)123 










CPPACh 


Wilson 


4.55(18) 


115(2) 


143(6) 


Mp 
\a — > \)) 


(lyyy jtij 










CP-FATS 
(2000)EI 


Iwasaki-l-SW 


4.4(2) 


110(4) 


132(6) 


Mp 
(a ^ 0) 










ALPHA^UKQCD 


0(a) SW 




97(4) 




Ik 


(1999) Ea 










(a ^ 0) 


QCDSE^ 


0(a) SW 


4.4(2) 


105(4) 




fo 


(1999)EI 










(a ^ 0) 


QCDSE, 


Wilson 


3.8(6) 


87(15) 




?'o 


(1999)Ea 










(a ^ 0) 


RIKENjBNL 


DWF 




105(6)(21) 


127(6)(25) 


Mp 


(2000)B 










a-^ « 2.1 


CP-PACS 
(2000)E9| 


n/ = 2 


o 44+0.14 


8811 




Mp 


Iwasaki+SW 








(a->0) 


QCDS£,-UKQCD 


n/ = 2 


3.5(2) 


90(5) 




fo 


(2000)y 


0(a) SW 








[1.9-2.2] 



masses, where the effects of the artifacts (enhanced chiral logs) becomes significant. 
So, in fact, the validity of quenched calculations will be judged a 'posteriori by 
comparison with (2 -I- 1) flavor calculations. 

Quenched LQCD simulations proceed with the assumption that the stable par- 
ticle masses and decay constants are affected at roughly the 10% level, i.e. the 
deviations are "small" . Within this assumption, the quenched results can be pre- 
sented in one of two ways: (Pi) analyze the data in all possible ways and take the 
variation in the estimates as a measure of the quenching uncertainty, or (P2) find 
combinations of quantities for which quenching effects are expected to be small, and 
use only these to extract a quenched number (which is then presumably "closest" 
to the real world). The results in Table || represent both approaches; for example, 
calculations by JLQCD and CF-PACS are closer to Vx, while ALFHA-UKQCD 
follow the Vi philosophy. 

fe start with three concensus statements regarding the quenched data in Ta- 
First, for a given gauge and fermion action, the raw lattice data, and the 
masses of hadrons as a function of the bare quark masses, are consistent between 
the different calculations. Second, there is ~ 10% uncertainty in the values of quark 
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masses depending on what physical quantity (Mp, Mk*, fx, /ir, or tq) is used to set 
the lattice scale. (In this case, for example, proponents of 7^2 would argue that Mp 
should not be used to set a because of its large width and similarly large shift in the 
mass between the unmixed and physical state.) Third, there is 20% uncertainty 
in nis depending on whether Mk or AIk* or is used to fix it {Mk* and M^ 
give consistent results). Furthermore, as expected,, neither ms{MK) nor ms{M^) 
yield the observed splittings in the baryon decuplet e3o. Our view, in summary, is 
that these variations, with the states used to fix a and the quark masses, provide a 
very useful handle for evaluating the "reality" of coming simulations with dynam- 
ical quarks: these differences should vanish as the input dynamical quark masses 
are tuned closer and closer to their physical values. 

Three other factors give rise to some of the differences observed in Table ||. These 
are (i) chiral extrapolations (ii) continuum extrapolations, (iii) the renormalization 
factor connecting lattice quark masses to those in the MS scheme at = 2 GeV. A 
brief discussion of these is as follows. 

Chiral Extrapolations: Simulations with physical masses for u and d quarks are 
computationally too expensive and beyond the capability of today's computers. 
Consequently one simulates the theory over a range of quark masses, typically 
2ms Ws/3, and extrapolates to to using predictions of (quenched, or partially 
quenched for nj = 2) chiral perturbation theory. 

Continuum extrapolations: The lattice theory at scale a has discretization errors 
of 0{ala'^), where i, n depends on the order of improvement of the lattice action 
and operators. These errors are removed by extrapolating results calculated over a 
range of scales, typically 2 < 1/a < 4 GeV, to a = 0. 

Renormalization Constants: These connect the lattice quark masses to either 
the renormalization group invariant mass lj, or those in a continuum renormal- 
ization scheme like MS at some fixed scale. These factors have been estimated 
using perturbation theory (CP-PACS), semi non-perturbatively (APE), and fully 
non-perturbatively (ALPHA-UKQCD); the latter being the most accurate. 

One way to improve the reliability of the chiral and continuum extrapolations 
is to use quantities (or their ratios) that have the least ambiguity and dependence 
on TO or a (7^2 approach). For example, while the statistical signal for extracting 
Mp is good, using it to set the lattice spacing can lead to a ~ 10% uncertainty 
since the mass of the stable p calculated in lattice simulations could differ by ^ 100 
MeV from the physical value. It is therefore desirable to use ratios with the bos4 
statistical signal, and the smallest variation with the quark masses or a (see Ref. 123 
for discussion). 

It is not straightforward to assess the absolute magnitude of these three errors, 
independent of the analyses. The reason is that they are correlated and depend on 
the details of how each analysis was done. For example, consider doing a continuum 
extrapolation with simulations at 1/a = 2, 3, and 4 GeV, scales typical of the 
best calculations. If the relative (systematic) error at the three points, due to say 
the chiral extrapolation, is not the same, then the bias can get magnified by the 
continuum extrapolation. From a survey of the various calculations, our estimate 
of the range of uncertainty in TOs from these three effects varies between 2 — 15 
MeV, with the lower value coming from the ALPHA-UKQCD result which includes 
a fully non-perturbative calculation of the renormalization constant, and considers 
a kaon composed of two quarks of mass « TOs/2 rather than to^ and to„. The 
latter simplification avoids the need for chiral extrapolation in to„ and is justified 
on the basis that the chiral fit, M^ versus to, shows a linear behavior in the range 
2tos — ms/2. 

Quenched calculations have significantly improved our understanding of, and 
ability to control, these three sources of systematic errors. For example, (i) use of 
improved actions has lead to a significant reduction in discretization errors (with 
even just 0{a) improvement, the residual error now is only 5 — 10% at 1/a ^ 2 GeV); 
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(ii) development of reliable methods for calculating the renormalizatitp. constants 
non-perturbatively has reduced the associated uncertaintly to < 2% iB; and (iii) 
partially quenched ChPT e2I has provided us with a much improved understanding 
of how to extract physical results from heavier quarks (in the range 2ms — "^8/4), 
especially in the case of future 2+1 dynamical flavor simulations. To summarize, 
our bottom line, based on quenched simulations, is 

• The technology for all aspects of the calculation has been developed. 

• We have a very good understanding of statistical and systematic errors. 

• The current best estimates are: 90 < < 140 MeV and 3.8 < m < 4.8 MeV. 

4. Lattice QCD: Results with 2 dynamical flavors 

Until very recently, the least controlled systematic uncertainty, as discussed 
above, was that due to the quenched approximation, which was forced on us by 
limitations in computer resources. As of this year a number of collaborations (CP- 
PACS, JLQCD, APE, QCDSF-UKQCD, MILC, SESAM) have begun to report 
results from simulations with two dynamical flavors. These were reviewed at Lattice 
2000 by Lubicz E3 and we discuss them briefly later. For the most part, we restrki 
our attention to the recent published results from the CP-PACS collaboration Ej 
(see Table as these results are of quality similar to the quenched calculations in 
terms of statistics, lattice sizes, and range of a values explored. 

The CP-PACS results are remarkable in that they show amazing consistency 
between ms{MK) and ms{M^), in contrast to the ~ 20% difference seen in the 
quenched theory. Also, compared to the quenched results of the ALPHA-UKQCD 
collaboration (who analyze data under the V2 philosophy), these are « 10% lower; 
the "expected" size of quenching errors. 

There, however, remain issues that need be be better understood. First, com- 
paring CP-PACS quenched and unquenched calculations with the same gauge and 
fermion action and with similar analysis methods, the unquenching effect is much 
> 10%. For example, ms{MK) changes from 110(4) 88tg MeV and ms{M^) 
from 132(6) 90tli MeV! Second, extrapolations in sea quark mass from "heavy" 
dynamical u and d quarks, i.e. in the range 2ms —>- yield already (within the 

statistical uncertainty) very good agreement between ms{Mx) and ms{M^). One 
needs to check, however, tbat this agreement persists when the partially quenched 
analyses recommended in e3 are fully implemented. Third, the role of the (so-far) 
neglected strange sea quarks remains to be investigated. Note that so far all lattice 
calculations indicate that the effect of adding dynamical flavors is to lower the es- 
timates of quark masses. Fourth, the renormalization constants used for all n/ = 2 
estimates are perturbative; the corresponding non-pcrturbative calculations have 
yet to be done. 

In spite of these cautionary remarks, we consider the CP-PACS numbers the 
current best lattice estimates. We quote these as [m„ + md](2 GeV) = 6.8±0.8 MeV 
and ms(2 GeV) = 89±11 MeV, where we have doubled the error on mu + m^ quoted 
by CP-PACS, based on that on rUg. 

Finally, a brief comment on how CP-PACS numbeiis-|Compare with preliminary 
data reported at LATTICE 2000 (seejieview by Lubicz E3) . The numbers vaxy from 
ms = 90(5) MeV (QCDSF-UKQCD El) to ms « 110 MeV (APE M, MILC 111). We 
expect that within the next couple of years a number of collaborations will have 
Uf — 2 flavor results of numerical quality similar to those by CP-PACS. We are 
therefore confident that a consistent picture of light quark masses will emerge soon. 
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